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Abstract 

The structure of a vortex in the inner crust of neutron stars is calculated within the framework 
of quantum mean field theory taking into account the interaction with the nuclei composing the 
Coulomb lattice. Making use of the results obtained with different nuclear interactions, the pinning 
energy, relevant in the study of glitches, is worked out. Quantal size and density dependent effects 
are found to be important. 
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I. INTRODUCTION 



According to standard theoretical models, the inner crust of a neutron star, where a 
lattice of nuclei is surrounded by a sea of free neutrons and relativistic electrons, lies in 
the range of densities between 4 xlO 11 and 1.6 xlO 14 g/cm 3 ( 0.0014 n and 0.57 n , 
where n is the nuclear saturation density, n = 0.16 fm -3 = 2.8 10 14 g/cm 3 ). The first 
microscopic study of this system was carried out by Negele and Vautherin. They determined, 
making use of a Skyrme-like functional of the energy, the isotopic composition as well as 
the distance between the (spherical) nuclei forming the Coulomb lattice Subsequent 
studies found that for densities larger than about no/3, other shapes rather than spherical 

n n 

become favoured [2]- [4]. This is still a debated topic. In any case, in our study we shall limit 
ourselves to densities n < no/3, where the assumption of spherical nuclei is considered to 
be safe. Another assumption at the basis of Negele and Vautherin's study was the Wigner- 
Seitz approximation, in which one studies a single cell of the lattice. Only very recently 
calculations using band theory have been performed j^]. 

Studies of neutron matter testify to the fact that the system is superfluid in a wide range 
of densities [6]. This result implies that the rotation of the star leads to the formation of 
vortices in the crust. They are expected to form an array whose average density per unit 
area is given by the Onsager-Feynman formula: 

P = 4m^, (1) 

where Q is the angular velocity of the star and m is the nucleon mass. The above relation 
leads to an average (macroscopic) distance between vortices of the order of 10~ 3 cm for 
Q ~ 1 ms _1 . 

The presence of impurities - represented, in the present case, by the nuclei constituting 
the lattice - has a profound influence on pairing correlations and on vortex dynamics. In 
fact, momentum is much higher inside the nuclear volume than in the surrounding neutron 
sea. This fact leads, as a rule, to a decrease of the pairing gap. Vortices, whose core is made 
out of normal matter, may thus find it energetically favourable to pin to the impurities, 
reducing the loss of pairing energy. 

One of the reasons for the interest in the interaction between vortex and nuclei in the 
inner crust is due to its relation with a possible explanation for the phenomenon of glitches, 
the sudden spinup of the rotation frequency observed in many neutron stars |l[ Q|. 



Considerable information concerning glitches has been accumulated over the years, and 
their properties are compatible with the inner crust as their origin J9]. However, a detailed 
comparison between observations and theoretical models has never been attempted. This 
is partly due to the fact that even the basic static interaction between a single nucleus and 
a vortex has never been studied in detail with state-of-the-art micoscopic theory, let alone 
the much more complex issue of vortex dynamics in the inner crust. Only semiclassical 
estimates of the pinning energy - the energy cost to build up a vortex on a nucleus, rather 

nnn 

than far from it - have been available [10, LUU, H2|], while the first quant al study of a vortex 
in uniform neutron matter was worked out only a few years ago [13l |. 

In the following we present a detailed account of a microscopic study, based on Hartree- 
Fock-Bogoliubov (HFB) quantum mean field theory, of the structure of a vortex in the inner 
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crust of a neutron star ( partial accounts of this work were published in refs. [U, H5|, H6|). In 
the present paper we aim at clarifying some basic qualitative features of the nucleus- vortex 
interaction, which, to our knowledge, have not been considered until now (in particular 
quantal finite size effects) and are absent from semiclassical models. We are well aware that 
our work represents only a starting point for a comprehensive understanding of vortices in 
the inner crust of neutron stars and of their possible relation to the phenomenon of glitches, 
because of various limitations which remain at the basis of our approach. In particular: 

- We treat pairing only at the mean field level. It is well known that medium polarization 
effects have a deep influence on pairing correlations in homogeneous matter [6||. This is even 
more so in the case of the inhomogeneous systems associated with each of the Wigner-Seitz 



cells studied by Negele and Vautherin 
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- We restrict ourselves to the region of spherical nuclei, while the largest contribution to 
the moment of inertia of the inner crust come from the deepest layers, where the assumption 
of spherical nuclei breaks down, and other phases (like the pasta phase) become important. 

- We consider the interaction between a single nucleus and the vortex, neglecting the 
influence of other nuclei building the Coulomb lattice. The soundness of this approximation 
depends on the vortex radius - becoming unacceptable when the vortex radius is of the 
order, or larger than, the radius of the Wigner-Seitz cell. 

- We assume axial symmetry for the nucleus-vortex system. As a consequence we can 
only consider special spatial configurations, and we are not able to calculate the pinning 
energy as a function of the distance between the center of the nucleus and the vortex axis. 
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- In our calculations we assume a simple isotopic composition for the nuclei, using always 
the same number of protons Z = 40, corresponding to a proton shell closure, as we neglect 
the spin-orbit term in the single-particle potential. A better calculation should consider 
the spin-orbit term of the mean field and use an updated version of the original Negele 
and Vautherin analysis in which the effects of pairing correlations are considered in 
determining the favoured number of protons and the lattice step, as a function of density. 

Finally, there are limitations which are intrinsic to the present status of the theory of 
nuclear interactions and nuclear structure, in particular concerning the choice of the ef- 
fective interaction used to produce the mean field. In fact, we shall see that some of our 
results depend, sometimes even qualitatively, on particular features of the selected nuclear 
interaction. 



II. HARTREE-FOCK-BOGOLIUBOV EQUATIONS 

We want to study the properties of nuclei embedded in a sea of superfluid free neutrons, 
as well as the properties of particular excitations of such a system, namely vortices. The best 
available tool to carry out such a study at the mean field level is provided by the Hartree- 
Fock-Bogoliubov (HFB) theory. The associated equations are widely used in both nuclear 
and solid state physics, where they are better known as the Bogoliubov-De Gennes equations 

HQ- 

The HFB equations allow, starting from two-body forces V p h and V pa i r (assumed to 
be of zero range) in the particle-hole channel and in the particle-particle channel respectively, 
to determine in a self-consistent way the single-particle Hartree-Fock Hamiltonian H = 
T + U HF and the pairing field A(x), as well as the amplitudes u i) Vi and the energies E { of 
the quasi particles. For a given particle species (neutrons or protons), the HFB equations 
can be written as 

(#(x) - AK(x) + A(x)^(x) = E iUi (x) 

A*(x)^(x) - (F(x) - AH(x) = £to(x) 
where A denotes the chemical potential, which is fixed so as to yield the desired average 
number of particles. 

While this framework is quite solid, the situation is very different concerning the choice 
of residual interactions acting in the particle-hole and in the particle-particle channel. The 
situation is made particularly trying, because in the Wigner-Seitz cell we must deal both 
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with strongly asymmetric matter (in the interior of the nuclei) and with low-density neutron 
matter. The parameters of the commonly adopted effective forces are usually derived fitting 
binding energies and other properties of atomic nuclei (mostly magic nuclei) in their ground 
state, as well as the binding energy of nuclear matter. In some cases also results in neutron 
matter have been taken into account, but at higher densities than those relevant here. On 
the other hand, many-body calculations should give reliable results for the equation of 
state of low-density neutron matter, because the uncertainties related to three-body force 
play a minor role. The density dependence of neutron matter energy resulting from such 
calculations has in fact been compared to that of current Skyrme forces (cf. e.g. j^l, Q)). 
We shall make use of the two Skyrme forces SkM* and SLy4, which display a number of 
attractive properties. The value of the energy per nucleon E/A, calculated with Bruckner 
theory based on a G-matrix derived by the bare nucleon-potential, and shown in Fig. El 
lies in between the results obtained with these two forces. A very relevant quantity in our 
studies appears to be the effective mass associated with the effective forces. To investigate 
this point, we have decided to show also results obtained with two other forces: SII and 
SGII, whose effective mass is close to that of SLy4 and of SkM* respectively (cf. Fig. CD). 
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Figure 1: (left) The energy per baryon calculated in neutron matter with Bruckner theory with the 
bare nucleon- nucleon Argonne interaction Vis (dot-dashed curve) [21] is compared to the results 
obtained with the Skyrme interactions SII (solid curve), SLy4 (dashed curve), SkM* (short dashed 
curve ) and SGII (dotted curve). 
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Figure 2: The effective mass in neutron matter (left panel) and in symmetric nuclear matter (right 
panel) for the different Skyrme interactions used in this paper: SII (solid curve), SLy4 (dashed 
curve), SkM* (short dashed curve), and SGII (dotted curve). 

Concerning the particle-particle channel, we shall use the density dependent contact force 



whose parameters V = -481 MeV fm 3 , 77 = 0.7 and a = 0.45 have been chosen so as to 
reproduce the pairi ng g ap of neutron matter obtained making use of a realistic nucleon- 
nucleon interaction [22||. Because of the zero-range nature of the interaction, all values of 
the momentum transferred in the nucleon-nucleon interaction processes are allowed. As 
a consequence one needs to introduce a cutoff parameter E cut . We have implemented this 
condition in the basis used to calculate the pairing matrix elements (cf. Appendix A), taking 
only single-particle states with energies lower than E cut = 60 MeV. In Fig. [3] we compare 
the pairing gaps at the Fermi energy calculated in neutron matter with this force and with 
the bare Argonne nucleon-nucleon potential. In Fig. [3](a) and (b) we use the single-particle 
levels calculated with the SLy4 and the SkM* interaction respectively: the latter yields 
higher values of the gap, because of the higher effective mass and thus of the higher single- 
particle level density. The contact force reproduces rather well the density-dependence of 
the gap obtained with the bare force, although it underestimates the value at the maximum 



The single-particle potential is obtained self-consistently from the calculated density 
n (x) = J2i v i( x -) 2 ' We shall follow ref. j^j, but we shall neglect the spin-orbit H so and 
spin coupling H sg terms in the Hamiltonian (cf. the comments at the end of Section IV B 




(3) 



by 10-20%. 
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Figure 3: (a) The BCS pairing gap Ajr in uniform neutron matter at the Fermi energy, is shown 
as a function of the Fermi momentum kp — (Stt 2 /!) 1 / 3 , calculated with the zero-range, density 
dependent pairing interaction ([3j) (dashed curve) or with the Argonne interaction (solid curve). In 
both cases the single-particle levels have been calculated with the SLy4 interaction, (b) The same, 
but with single-particle levels calculated with the SkM* interaction. 



nd Fig. 1211) . Furthermore, we shall neglect the terms associated with the time-odd momen- 
turn density (see for example [26]), which are associated with the presence of the vortex. In 
fact, we have evaluated their contribution to the energy perturbatively, finding that their 
contribution is negligible for the scope of this paper. 

The pairing gap is obtained self-consistently from the so-called abnormal density ft(x) = 
\ Y^i ^( x K*( x )> according to 

A(x) = -V^ r (x)/<x). (4) 

The HFB equations are solved enclosing the system in a cylindrical box of radius pb ox 
= 30 fm and height hb ox = 40 fm (cf. Fig. [4]). We assume perfectly reflecting walls, with 
single-particle wavefunctions which vanish at the edges of the box (see the remarks on this 
issue in Appendix A). 

In the calculations discussed in the following, we shall only consider axial symmetry, in 
which case the pairing field can be written as jsifl, Q, 3) 

A(p,z,<P) = A(p,z)e^, (5) 

v being the vortex number, which indicates the number of quanta of angular momentum 
carried by each Cooper pair. In particular, if we set v = we obtain the usual HFB 
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Figure 4: Geometrical features of the cylindrical cell used in most calculations of this paper. 

equations (no vortex). Setting instead v — 1, we can describe a vortex excitation, in which 
each Cooper pair carries one unit of angular momentum along the z— axis. Experiments 
on superfluids indicate that it is energetically more favorable to develop an array of v — 1 
vortices rather than a few vortices carrying many quanta of angular momentum and in the 
following we shall apply this condition also in neutron stars. When v ^ 0, one obtains a 
current given by 

J{r,z,j>) = --—l^ v i^ z ^) al ' 6 

mp ~ J ocp 

and an associated velocity field J /in. Aside from axial symmetry, the system has mirror 
symmetry respect to the x — y plane (see Fig. BJ, so that A(p, z) = A(p, —z). Applying the 
parity operator (0 — > (j) + tt and z —z) to the pairing field we obtain 

PA(p, z)e iv * = e^A(p, z)e^. (7) 

Therefore the pairing field is an eigenstate of the parity operator with eigenvalue (— l) u . 
This condition implies that the Cooper pairs involved in the v — 1 excitation are constructed 
out of single particle levels of opposite parity. 
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Rather than solving the HFB equations directly in coordinate space, we expand them on 
a single-particle basis. For details concerning the numerical procedure, we refer to Appendix 
A. 

III. THE PINNING ENERGY 

Models that relate the glitch phenomenon to the vortex dynamics in the inner crust, 
usually assume that the vortex lines which thread the neutron superfluid pin to the crust, 
being attracted or repelled by the nuclei forming the Coulomb lattice. If the vortex lines 
remain fixed, the velocity of the superfluid remains constant, while the velocity of the crust 
decreases due to the magnetic braking. At some critical velocity difference, the Magnus 
force unpins the vortex lines, and angular momentum is given back to the star, causing the 
glitch. 

The determination of an 'optimal' vortex line depends, besides the lattice properties and 
the vortex tension, on the vortex-nucleus interaction, or, more precisely, on the free energy 
as a function of the distance between the nucleus and the vortex axis. As it was mentioned 
in the Introduction, although in the present paper we shall improve in different ways on 
previous models, we shall only be able to compare the two limiting configurations: a vortex 
whose axis passes through the center of the nucleus, and a vortex placed far from the nucleus. 

We shall now discuss more precisely how we compute the pinning energy, E pinningi at a 

Bm density in the crust, following (except for a sign) the definition by Epstein and Baym 
. At a given density, they calculated, within the Ginzburg-Landau approximation, the 
energy of the vortex-nucleus system as a function of the distance between the vortex axis 
and the center of the nucleus, determining the configuration of minimum energy, which 
often turned out to correspond to zero distance. Then they calculated the pinning energy 
as the difference between this minimum energy and the energy of a configuration where 
the vortex and the nucleus are sufficiently far apart, so that their mutual interaction can 
be neglected. Due to the axial symmetry condition, we shall only compare the energies of 
the zero-distance configuration (vortex axis passing through the center of the nucleus) and 
of the non-interacting one. In this framework we only deal with one nucleus at a time, 
neglecting the interaction of the vortex with distant nuclei. The situation can then be 
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Figure 5: (left) The pinning energy E p i nning is obtained taking the difference between the energy of 
the pinned configuration (1), in which the vortex axis passes through the center of the nucleus, and 
the energy of the interstitial configuration (2), in which the vortex is far from the nucleus, (right) 
We rearrange the configurations on the left, indicating that the pinning energy is equivalent to the 
difference between the energy cost to build the vortex on a nucleus (E^ c ) and in uniform matter 

( E uZif)> S0 tnat E pinning = _ E unif 

schematically depicted as in the left part of Fig. [51 We have to compare the energies of 
the pinned configuration (1) (vortex axis passing through the center of the nucleus) and of 
the interstitial configuration (2) (vortex axis far from the nucleus). We can then consider 



pinning 



as the difference between two excitation energies: 



pinning 



rpvor rpvor 

nuc unif > 



(8) 



where E^j is the energy of a vortex in the uniform matter configuration, compared to the 
energy of the same configuration in its ground state, while E^[ c is the energy of the vortex 
pinned on a nucleus (immersed in the neutron sea), compared to the energy of the isolated 
nucleus immersed in the neutron sea; that is, 



rpvor 771 771 

nuc -^pinned -^nucleus 

rpvor 77 77 

unif ^vortex -^uniform) 



(9) 



where the energy associated with each of the four configurations is calculated as described 
in Appendix A. Pinning the vortex on the nucleus is energetically convenient, if the pinning 
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energy is negative (note that Epstein and Baym define the pinning energy with the opposite 
sign). We note that the total and pairing energy of the various configurations, as well as 
the energy cost to create the vortex, depend on the size of the box. On the contrary, the 
value of the pinning energy, which is the result of a phenomenon localized on the nucleus, 
converges for a sufficiently large box. 

In practice, four different HFB calculations have to be performed, using the same cylin- 
drical cell (cf. Fig. [6]): one with Z = 40, v = 1 (Fig. E](a)), yielding E pinnec i; one with 
Z = 40, v = (Fig. 0(b)), yielding E nudeus ; one with Z = 0, v = 1 (Fig. 0(c)), yielding 
Evortex] and finally, one with Z = 0, v = (Fig. E](d)), yielding E uni f orm . To obtain the 
correct pinning energy it is essential that the calculations refer to the same asymptotic neu- 
tron density. In the calculations (b) and (d), without the vortex, we use the same value of 
the chemical potential A, yielding N nuc i eus and N uni f orm neutrons in the cylindrical cell. The 
presence of the vortex in calculations (a) and (c) leads in each case to a small decrease of 
the number of particles. We compensate this reduction by a slight increase in the value of 
A. In practice, rather than attempting a very fine tuning of A we prefer to account for the 
residual difference in the number of particles subtracting the terms 

^E nuc A ( ^'pinned ^nucleus) > ^E un j j j A (A/ T vortex ^uniforrn) (^^) 

respectively to E^ c and to E^j. Even if the pinning energies represent only a small fraction 
of the total energy, of the order of 10 -3 — 10 -4 , the subtraction scheme we have just outlined 
produces numerically reliable results (cf. Appendix B), which will be presented in the next 
Section. One should notice that the size of the cylindrical cell does not have to coincide 
with that of the physical Wigner-Seitz cell; it must only be large enough, so as to obtain 
convergent results for E pinning , It is clear, however, that neglecting neighbouring nuclei can 
be inconsistent, if the radius of the box becomes of the order of the lattice constant. This 
point will be further discussed below in Section IV C. 



IV. RESULTS 



We have performed calculations at different densities in the inner crust, ranging between 
n « 0.001 fm -3 and n ~ 0.04 fm -3 . Our discussion will be mostly based on the results 
obtained at eight densities with the interaction SLy4 and SkM*, but we shall also present 



11 




Figure 6: Schematic picture of the four different physical configurations needed to calculate the 
pinning energy: (a) cell with a vortex pinned on a nucleus ; (b) cell with a nucleus in the neutron 
sea ; (c) cell with a vortex in the neutron sea; (d) cell in the uniform neutron sea. 

results obtained with the interactions SII and SGII (cf. Tables I and II). 

In the following we first present the results obtained for a nucleus embedded in the neutron 
sea, in the absence of vortex (HFB calculations with v — 0). We then discuss the structure 
of the vortex [y — 1), without or with the nucleus. Finally, we present the pinning energies 
obtained as a function of density. 
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Interaction 


«oo [fm' 3 ] 


a [fm] 


n [fm~ 3 ] 


R [fm] 


SLy4 


0.001 


0.82 


0.098 


6.06 


SLy4 


0.002 


0.82 


0.098 


6.04 


SLy4 


0.004 


0.82 


0.090 


6.81 


SLy4 


0.008 


0.86 


0.091 


6.75 


SLy4 


0.011 


0.91 


0.089 


6.76 


SLy4 


0.017 


0.96 


0.081 


7.33 


SLy4 


0.026 


1.05 


0.070 


7.83 


SLy4 


0.037 


1.13 


0.057 


8.43 


SII 


0.001 


0.73 


0.093 


6.15 


SII 


0.002 


0.73 


0.093 


6.14 


SII 


0.004 


0.74 


0.092 


6.10 


SII 


0.008 


0.77 


0.086 


6.74 


SII 


0.011 


0.78 


0.085 


6.71 


SII 


0.016 


0.81 


0.082 


6.75 


SII 


0.024 


0.87 


0.073 


7.18 


SII 


0.035 


0.94 


0.062 


7.60 



Table I: Parameters of the Fermi function (fTTj) . fitting the calculated neutron densities in the Wigner 
Seitz cells with the SLy4 and SII interactions. The values of the density have been rounded off at 
the third digit. 
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Interaction 


«oo [fm~ 3 ] 


a [fm] 


n [fm~ 3 ] 


R [fm] 


SkM* 


0.001 


0.88 


0.097 


6.12 


SkM* 


0.002 


0.90 


0.097 


6.17 


SkM* 


0.004 


0.96 


0.091 


6.76 


SkM* 


0.008 


0.97 


0.088 


6.73 


SkM* 


0.012 


1.00 


0.086 


6.78 


SkM* 


0.016 


1.05 


0.080 


7.03 


SkM* 


0.025 


1.09 


0.070 


7.33 


SkM* 


0.038 


1.17 


0.054 


7.95 


SGII 


0.001 


0.86 


0.097 


6.91 


SGII 


0.002 


0.88 


0.097 


6.91 


SGII 


0.004 


0.90 


0.096 


6.82 


SGII 


0.008 


0.94 


0.092 


6.86 


SGII 


0.011 


1.00 


0.086 


7.34 


SGII 


0.016 


1.00 


0.082 


7.38 


SGII 


0.026 


1.07 


0.070 


7.77 


SGII 


0.037 


1.13 


0.057 


8.27 



Table II: Parameters of the Fermi function (fTTj) . fitting the calculated neutron densities in the 
Wigner Seitz cells with the SkM* and SGII interactions. The values of the density have been 
rounded off at the third digit. 
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A. Calculations with v = 



We first discuss the density and the pairing gaps obtained solving the HFB equations 
for v = in the cylindrical box, for the nuclei immersed in the neutron sea at the various 
densities. The pairing field in the Wigner-Seitz cell has been studied by various groups, 
either based on Woods-Saxon potentials {27] or on selfconsistent HFB calculations with 



Skyrme forces [18, 28, |29j. The effects of temperature, and the specific heat, have also 
been discussed [27|, |29j, but we shall limit ourselves to zero temperature. Our calculations 
are performed assuming axial symmetry, and they should respect the spherical symmetry of 
isolated nuclei. In our calculations the neutron density is a function of the coordinates p and 
z. To test the sphericity of our calculations we evaluate the mean value of the density within 
a spherical shell 0.05 fm wide and centered around a mean distance r from the center of the 
nucleus and then evaluate the fluctuations at every point. The relative error (rms/average 
density) is everywhere lower than 0.03. 

As first remarked by Negele and Vautherin [1], nuclei in the inner crust turn out to be 
fatter than atomic nuclei, due to their interaction with the neutron sea. We have fitted the 
neutron density resulting from our calculations with a Fermi function: 

/W = "~ + 7 \ (H) 



(i + «p(t=ai)) 



where gives the asymptotic value of the neutron density far from the nucleus, and we 
have determined the values of n , Ro and a which best reproduce the computed density 
profile and are reported in Tables I and II. In the calculations we have adjusted the value of 
A so as to obtain the same set of values of for the various interactions, but there remain 
slight differences in a few values (cf. also Tables III and IV). We compare the calculated 
and fitted densities for three cases in Fig. [3 The diffusivity a increases from about 0.7 to 
about 1.1 fm with increasing density, to be compared with the typical value of 0.65 fm for 
finite nuclei, and the radius R (up to 8 fm at the highest densities) also turns out to be 
rather large compared to ordinary nuclei. 

In Figs. IHliTOl we show the pairing gaps and selfconsistent potentials as a function of the 
distance from the nucleus, calculated for w 0.001, 0.01 and 0.035 fm -3 ; the associated 
effective masses are shown in Fig. dU 

As we anticipated (cf. FigEj), the behaviour of the forces can be divided into two groups: 

15 



0.12 

0.1 y&** 
0.08 
0.06 
0.04 
0.02 




V 
V 



i^-mfH+iiii 



0.12 



10 12 14 



p[fm] 



0.1 
0.09 
0.08 
0.07 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 





10 12 14 



p[fm] 



Figure 7: The neutron density calculated in the cylindrical box with the nucleus on the z — plane 
(crosses) is compared with the fit obtained with the Fermi function (cf. Eq. [TTj) . As we discuss 
in the text, the calculated density is spherically symmetric to a good approximation. The figures 
refer to the asymptotic densities noo= 0.001, 0.011 and 0.037 fm -3 , and have been computed with 
the SLy4 interaction. 

the effective masses associated with SGII and SkM* are quite similar ( rrf jm w 0.9 in the 
interior of the nucleus), as are those associated with SLy4 and SII {rrf jm w 0.7). In order 
to reproduce the same asymptotic density the SII and SLy4 interactions produce a deeper 
selfconsistent potential than the SkM* and SGII. In all cases the pairing gap increases, going 
from the nuclear region to the sea of superfluid neutrons. This is in keeping with the density 
dependence of the pairing gap in uniform matter, previously shown in Fig. 02 the interior 
of the nucleus corresponds to large values of the local Fermi momentum where the pairing 
gap tends to vanish, while even at the smallest densities considered here, the gap in neutron 
matter is at least close to 1 MeV. In the case of the SkM* and SGII interactions the pairing 
gaps show a peak in the nuclear surface region for w 0.001 fm -3 (cf. Fig. Ej). The pairing 
gaps are smaller with SII and SLy4, due to their lower effective mass and level density. 
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Figure 8: Pairing gap (left) and selfconsistent potential (right) in the cylindrical cell without the 
vortex [Z = 40, v = 0) for « 0.001 fm -3 (hp ~ 0.3 fm -1 ). The different interactions are: SLy4 
(dashed curve), SkM* (short dashed curve), SII (solid curve), and SGII (dotted curve). 




p [fm] P [ fm l 

Figure 9: Pairing gap (left) and selfconsistent potential (right) in the cylindrical cell [Z = 40, v = 0) 
for « 0.011 fm -3 (hp ~ 0.7 fm -1 ). The different interactions are: SLy4 (dashed curve), SkM* 
(short dashed curve), SII (solid curve) and SGII (dotted curve). 




p [fm] P [ fm l 

Figure 10: Pairing gap (left) and selfconsistent potential (right) in the cylindrical cell (Z = 40, v = 
0) for rioc « 0.035 fm -3 (hp « 1 fm -1 ). The different interactions are: SLy4 (dashed curve), SkM* 
(short dashed curve), SII (solid curve), and SGII (dotted curve). 
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Figure 11: Effective masses associated with the calculations presented in Fig. [EK^oo ~ 0.001 fm -3 , 
top left), in Fig. [9] (n^ « 0.01 fm~ 3 , top right), and in Fig. [TO] (n^ « 0.035 fm~ 3 , bottom) for 
the four interactions SII (solid curve), SLy4 (dashed curve), SkM* (short dashed curve) and SGII 
(dotted curve). 

B. Calculations with v = 1 



We now turn our attention to the HFB calculations of the vortex [y = 1). We first 
discuss the results we obtain for the uniform neutron sea (Z — 0). In this case the vortex 
configuration keeps cylindrical symmetry, and a cut at a fixed value of z provides complete 
information. The density profiles and the pairing gap obtained with the SLy4 interaction 
at four different densities are reported in Fig. [121 The pairing gap vanishes on the vortex 
axis, and grows linearly for about 3-5 fm; then it slowly increases towards the value in v = 
neutron matter. In keeping with Fig. [3l the asymptotic value of A reaches its maximum 
for kp ~ 0.8 fm -1 , that is, for n ~ 0.02 fm -3 . Near the vortex axis the pairing gap displays 
oscillations that have been interpreted as Friedel-like oscillations, due to the presence of 
bound states in the vortex core J30I]. Our results are very similar to those obtained by 
Bulgac and Yu, who have solved the HFB equations for a vortex in uniform, infinite neutron 
matter for the first time [13]. In particular, the depletion of the density close to the vortex 
axis (cf. Fig. d2j) and the behaviour of the velocity field (cf. Fig. [13]) are in good agreement 
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with their results. Note however that we use a different pairing interaction and that we 
consider densities larger than those calculated in ref. [131 ]. 
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Figure 12: (left) Neutron density profiles for the vortex solution, without the nucleus at the center 
of the cell [Z = 0, v = 1), calculated with the SLy4 interaction at different densities: = 
0.0013 fin" 3 ; (solid curve, the density shown has been multiplied by 5), =0.011 fm 3 ; (dashed 
curve), rtoQ — 0.026 fm -3 (short dashed curve) and — 0.037 fm -3 ; (dotted curve). (Right) The 
corresponding pairing gaps. 
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Figure 13: Velocity field of a vortex calculated in the cell without the nucleus [Z = 0, v = 1), with 
the SLy4 interaction calculated with = 0.0013 fm ~ 3 ; (solid curve), = 0.011 fm ~ 3 ; (dashed 
curve , 0.026 fin -3 ; (short dashed curve , and 0.037 fm~ 3 ; (dotted curve) . For large values of p the 
velocity tends to the Onsager limit h/2mp. 



19 






Figure 14: The level curves (values in MeV) of the pairing gap associated with a vortex pinned on 
a nucleus [Z — 40, v = 1, top left), with a vortex in uniform matter [Z = 0, v = 1, top right) and 
with a nucleus in the neutron sea (Z = 40, v — 0, bottom) obtained with a SLy4 force for = 
0.011 fm ~ 3 . The nucleus is at p — z — and the vortex is along the z-axis. 

We now discuss the calculation of the vortex pinned on the nucleus (Z = 40, v = 1). The 
main results are illustrated in Figs. [H] and [151 which refer to a calculation performed at 
rioQ = 0.011 fm -3 with the SLy4 interaction. In Fig. [14] we show the gap in the p, z plane. 
In Fig. [15] the gap is shown as a function of p in the z = equatorial plane. In Fig. [15] we 
also make a comparison with the pairing gaps obtained for a vortex in uniform matter and 
for a nucleus in the neutron sea, previously discussed. The neutron density and the velocity 
field are shown in Fig. [161 

For large values of z the properties of the vortex are the same as those calculated for uni- 
form matter (cf. Fig. [14]). Approaching the nuclear surface the vortex is strongly modified, 
and the associated pairing gap is suppressed within its volume and also on its surface (for 
example A ^ 0.3 MeV at p w 7 fm and z =0 ), both compared to the case of an isolated 
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Figure 15: (left) The pairing gap associated with the pinned vortex on the z = plane [Z = 
40, v = 1, solid line) is compared with the gap associated with the nucleus in the absence of the 
vortex (Z = 40, v = 0, dashed line), (right) The pairing gap associated with the pinned vortex 
(Z = 40, v = 1, solid line) is now compared to the gap associated with the vortex in uniform matter 
[Z = 0, v = 1, short dashed line); the latter is also shown translated by 7.8 fm (dashed line). The 
calculation is for the SLy4 interaction at = 0.011 fm~ 6 . 
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Figure 16: (left) Density with (Z = 40, v = 1, solid line) and without (Z = 40, v = 0, dashed line) 
the vortex, calculated at = 0.011 fm -3 on the z = plane with the SLy4 force, (right) Velocity 
field of the pinned vortex. The different curves refer to different cuts; the dot-dashed line is the 
h/2mp function. 
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nucleus in the absence of the vortex (A « 1 MeV) and compared to the case of a vortex in 
uniform neutron matter (A « 2 MeV). The typical linear rise of the gap away from the vortex 
axis is delayed by about 8 fm, which is about the radius of the nucleus, as compared to the 
uniform case (cf. Fig. [15]). In other words, the vortex opens up and surrounds the nuclear 
volume, being unable to penetrate into it. The suppression of the v — 1 pairing field inside 
the nucleus was already found (but not thoroughly discussed) in a quantum calculation by 
Elgar0y and De Blasio jso|, although their study was not self-consistent and assumed an 
unrealistic cylindrical shape for the nucleus. The same effect appears concerning the velocity 
of the vortex flow. The velocity vanishes inside the nuclear volume, and the typical rise of 
the velocity field from zero to the asymptotic Onsager dependence v = h/2mp is delayed, 
again by about 8 fm at z = 0. The vortex induces some changes on the nuclear density, 
which becomes somewhat elongated along the vortex axis, increases close to the center of 
the nucleus and decreases close to the surface, mostly around the equatorial plane (cf. Fig. 
[T6l) . This behaviour can be compared with the case of the vortex in uniform matter, where 
the depletion takes place close to the axis (cf. Fig. fT2l above). 

The pairing gaps calculated for the four configurations at ~ 0.001, 0.01 and 0.035 
fm -3 with the interactions SLy4 and SkM* are collected in Figs. fT7tfT9l The results obtained 
with SLy4 at all densities reflect the main features discussed above (the results obtained with 
the SII interaction are very similar). One can observe that at the highest density, the results 
obtained with the SkM* interaction show a different behaviour. In fact, in this case the 
pairing gap associated with the vortex is not suppressed in the interior of the nucleus, but 
it shows a constant increase: the vortex can penetrate the nucleus. Differences are present 
also in the density (which presents a slight depletion close to the vortex axis) and in the 
velocity (which does not go to zero in the interior of the nucleus), which are closer to the 
case of the vortex in uniform matter. The pairing gap associated with the pinned vortex 
obtained with the SkM* (and with the SGII) interaction can be approximated rather well 
assuming that the gap suppressions produced by the nuclear field (A nuc i eus (p) / A^) and by 
the vortex in uniform matter (A vortex (p) / A^) respect to the value in uniform matter A^ 
act simultaneously and independently (cf. Fig. [20]). The superposition of the effects is then 
the product of the two suppressions: 

A(p) = (A vortex (p)/A OQ ) • (A nud (p)/A OQ ) • Aoo. (12) 
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This approximation would fail badly in the case of SLy4 interaction, where the pairing of 
the pinned vortex is suppressed much more strongly in the interior, as discussed above. 





Figure 17: The pairing gaps associated with the four different configurations (pinned vortex [Z — 
40, v = 1): solid curve; nucleus in the neutron sea (Z = 40, v = 0): dashed curve; vortex the 
neutron sea (Z = 0, v = 1): short dashed curve; uniform neutron sea [Z = 0, v = 0: dotted curve) 
calculated with the SLy4 (left) and the SkM* (right) interactions at ~ 0.001 fm -3 . 
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Figure 20: The pairing gaps calculated with the SGII (left) and the SkM* (right) interactions at 
rioo ~ 0.035 fm -3 (solid lines) are compared with the expression (jT2j) (dashed lines). 

The different behaviours discussed above can be understood qualitatively in terms of 
nuclear shell effects. As we remarked previously, the v — 1 Cooper pairs are formed out of 
single-particle levels of opposite parity. Shell effects strongly hinder the formation of such 
Cooper pairs inside the nuclear volume, as they separate the levels of opposite parity lying 
within the nuclear potential (i.e. the resonant states) by several MeV, even if these levels lie 
in the continuum. This can be seen in Fig. [221(a) and (b), where we show a calculation of the 
phase shifts 8\[E\ for / = — 10, in the mean field (without the vortex) associated with the 
interaction SLy4 at the densities = 0.011 fm~ 3 (cf. Fig. [9]) and = 0.037 fm~ 3 (cf. Fig. 
[TOll . In the figures the energy E is referred to the asymptotic value of the self-consistent 
potential U HF far from the nucleus (cf. Fig. [9] and Fig. [10]). In Fig. [22](c) and (d) we 
show the same calculation for the SkM* interaction. We recall that the resonant states and 
the resonant energy are characterized by the conditions 8i(E res ) = (2n + l)7r/2,n = 0, 1, .., 
and d8i(E res ) / dE > 0, which in the present cases are fulfilled by states with energies up to 
about 30-40 MeV. At low density (cf. Fig. [22](a) and (c)) the phase shifts are similar for 
the two interactions. Close to the Fermi energy A — V^, which is about 10 MeV, one finds a 
rather narrow resonance for I = 7 (T = tt / \d8i(E res ) /dE) w 3MeV ); one also finds narrow 
resonances for I = 6 (T w 0.5MeV ) and rather well defined resonances for I = 8 (T w 15 
MeV ), but they lie about 7 MeV below or above the Fermi energy. 

In Fig. [221(b) and (d) we show a calculation at higher density, corresponding to A — 
Vqo about 25 MeV. One still finds resonant states in the region around the Fermi energy, 
associated with / = 8 — 9, but they are much broader for the SkM* (the width becoming 
larger than the centroid energy) than for the SLy4 interaction (T w 9 MeV compared to 
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Figure 21: Schematic figure showing that the spin-orbit interaction tends to shift the energy of 
single-particle pairs involved in the formation of S = 0, v = 1 Cooper pairs by the same amount 
(cf. text). 

E res — Voq « 15 MeV for I = 8, and Y is still of the order of E res — for / = 9), in keeping 
with the fact that the mean-field potential is much shallower. The role of shell effects in the 
SkM* case is correspondingly much reduced, and the pairing interaction is able to mix levels 
of opposite parity and to create a pairing field of v — 1 character inside the nuclear volume. 

This interpretation is reinforced by a calculation of a v = 2 vortex. In this case, Cooper 
pairs are formed out of single-particle states of the same parity, and one can see in Fig. [23] 
that already for = 0.011 fm -3 the pairing gap obtained with the SLy4 interaction does 
not vanish within the nucleus: its value can be reproduced by the simple approximation of 
Eq. ( IT2l) . similar to the case of the SkM* interaction at high density (cf. Fig. 1201) . 

It could be argued that the spin-orbit interaction, not considered in the present calcu- 
lation, will shift the relative energies of levels of different parity around the Fermi energy, 
eventually modifying the pairing gap we obtain. However the v = 1 pairing field con- 
nects single-particle levels with opposite spin projections, and with L z -values such that 
L z (l) + L z (2) = v. This implies (except for L z = — 1) that the products L z (l) * S z (l) and 
L z (2) * S z (2) have the same sign, so that the two levels will be shifted by about the same 
amount by the spin-orbit interaction, as illustrated in Fig. EU One thus does not expect a 
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Figure 22: Phase shifts associated with the Hartree-Fock potential for a spherical nucleus, obtained 
with the SLy4 interaction for n^ = 0.011 fm -3 ( corresponding to A = 5.8 MeV (a)) and for n^ = 
0.037 fm -3 ( corresponding to A = 11.3 MeV (b)), and with the SkM* interaction for n^ = 0.011 
fm -3 ( corresponding to A = 5.8 MeV (c)), and for n^ = 0.038 fm -3 ( corresponding to A = 11.3 
MeV (d)). The energy is referred to the asymptotic value of the mean field potential far from the 
nucleus. 

significant change of the phase space density relevant for the vortex formation. 
C. Pinning energies 

We now turn to the calculation of the pinning energy, defined above in Section III. We 
recall that we compare the energy of the 'pinned' configuration (vortex axis passing through 
the center of the nucleus) with the 'interstitial' configuration, in which we assume that 
the vortex is placed far from the nuclei, so that its properties are the same as in uniform 
matter. In both cases, we assume that the interaction of the vortex with 'distant' nuclei can 
be neglected. In order to address the validity of our approximation in a quantitative way, 
one would need to perform systematic calculations taking into account the band structure 
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Figure 23: (left): The pairing gap associated with a v = 2 vortex pinned on the nucleus (Z = 
40, v = 1), calculated with the SLy4 interaction at = 0.011 fm -3 . The cuts refer to fixed value 
of z = 10 fm (short dashed curve), z = 5 fm (dashed line) and z = fm (solid line), (right) The 
pairing gap at z = (dashed line) is compared to the approximation (fT2|) (solid line) . 

associated with the lattice One expects, however, that our approximation should be 
reasonable if the vortex radial dimension is smaller than the distance between nuclei, as 
given by the diameter of the Wigner-Seitz cell calculated by Negele and Vautherin. To 
obtain some insight about this point, we need an estimate of the extension of the vortex. 
Different definitions of the vortex radius in uniform matter have been introduced. One often 
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Figure 24: Comparison between different definitions of the vortex radius: i?5o% (solid curve), i?go% 
(dashed curve) and coherence length (short dashed curve), calculated with the SLy4 interaction for 
a vortex in uniform neutron matter at various neutron densities. 



employs the coherence length, £ = f^kjp/mTvA^ shown in Fi g. I24l by the short-dashed curve. 
We also show two different estimates of the vortex radius [30||, namely the distances Rgo% 
or i?5o% from the axis where the pairing gap reaches 90% or 50% of its asymptotic value. 
According to this figure, the coherence length and Rqq% display a similar dependence on 
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density, while the quantity i?5o% is small and almost independent of the density. We shall 
take i?go% as a measure of the vortex radius, as this is the most conservative estimate. 
Let us then consider the interstitial configuration first. 

In this case we require that when the vortex has reached its asymptotic gap, the neutron 
density is still close to its asymptotic n^/^, that is, it is not appreciably affected by the 
neighbouring nuclei. We then introduce the distance R sa fe = Rws — 1.2Rn- According to 
the parametrization ( fTTl) . the quantity R sa fe represents the distance where the difference 
between the neutron density and its asymptotic value has reduced to 10% of its maximum 
value. It turns out (see Fig. [26]) that the dimension of the vortex is always smaller than 
Rsafe-) although it approaches this value at the largest densities we have considered in our 
calculations. 

Concerning the pinned configuration, the radii of the pinned vortices, measured through 
i?90%, are 4-6 fm larger than in uniform matter, depending on the interaction (cf. Fig. 1251) . 
However, the neglect of neighbouring nuclei appears to be even better justified for the pinned 
than for the interstitial configuration discussed previously. In fact, in the pinned case one 
should compare Rqo% with the distance between neighbouring nuclei 2Rws (or better with 
2R sa f e ). In all cases the calculated radii of the pinned vortices are considerably smaller than 
this value. 

Consequently we can safely describe the structure of the vortex in the inner crust, both 
in the interstitial and in the pinned configuration, within cylindrical boxes of the order of 
Rqo%. However, in the calculation of the pinning energy we need a common box size value 
for the two vortex configurations. This implies using a box of radius close to the value Rqo% 
for the pinned configuration, which is the largest one: this quantity remains smaller than 
R sa fei as can be seen comparing Fig. [25] with Fig. [261 

In practice we have used somewhat larger radii in the calculation of the pinning energy in 
order to check that our boundary conditions do not produce any relevant effect (see appendix 
A). 

The fact that the presence of neighbouring nuclei is not relevant for the calculation of the 
pinning energy for the densities considered in this paper, is confirmed by explicit estimates 
carried out in refs. (see Appendix C). 



28 




0.01 0.02 0.03 0.04 0.01 0.02 0.03 0.04 

rioo [fm" 3 ] n M [fm" 3 ] 



Figure 25: The vortex radius, Rqo% is shown for the pinned vortex (solid curves) and for the vortex 
in uniform matter (dashed curves), for the four Skyrme interactions. 

After calculating the energies of the four different configurations of Fig. El and using 
Eqs.([8j) and ([9j) we compute the pinning energies shown in Fig. EH for different values of 
A and (cf. Tables I and II). An estimate of the numerical errors associated with the 
calculations is given in Appendix B. Some details of the energy calculations performed with 
the interactions SLy4 and SkM* are reported in Tables III and IV, in which we list the 
values of the chemical potential, of the number of neutrons, of the total energy E tot and 
of the pairing energy E pair (cf. Appendix A). The total energy associated with the two 
calculations with a vortex (labeled 'Pinned' and 'Vortex'), has been corrected by the term 
AE introduced at the end of Section III, in order to account for the difference in number of 
neutrons compared to the corresponding calculations without vortex (labeled 'Nucleus' and 
'Uniform'). 

All the interactions yield weak pinning (E pinning < 0) at the smallest densities we have 
calculated (n^ ~ 0.001 fm -3 ), and yield a few MeV antipinning for = 0.01 — 0.02 fm -3 . 
At higher density, the SkM* and SGII interaction lead to a strong antipinning (E pinning > 
0) , in marked contrast with the other two interactions, which instead again produce weak 
pinning at the highest density points. This different density dependence can be related to 



29 




0.01 0.02 0.03 

n M [fm~ 3 ] 



0.04 



0.12 
0.1 



0.04 



I I I 








- R " \ 




100% 


- L2R » \ 






10% 






1 1 1 



10 
P [fm] 



15 



20 



Figure 26: (left) Comparison between the vortex radius (i? 90 %) and the distance R sa fe- The 
interactions SII and SLy4 have very similar radii at a fixed neutron density and their curves are 
almost superimposed; the same is true for SkM* and SGII. (right) Illustration of the density profile 
obtained with the parametrization (fTTj) . At a distance 1.2 Rq from the center of the nucleus, the 
difference between the neutron density and its asymptotic value has reduced to 10% of its maximum 
value. 

the qualitative differences in the vortex structure, discussed in the previous section. 
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Figure 27: ((left) Pinning energy calculated as a function of the asymptotic neutron density with 
the interactions SLy4 (dashed line) and SII (solid line). The lines connect the calculated points, 
(right) The same for the interactions SkM* (solid line) and SGII (dashed line). 
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Table III: Calculation of the pinning energy with the interaction SLy4 for the various neutron 
densities we have considered, (cf. Table I). In each case, we give the values of the chemical potential 
A, (in MeV), the number of neutrons N, the pairing energy E pa i r (in MeV), the total energy E tot 
(in MeV) (cf. Appendix A), associated with the four configurations labeled Pinned (vortex on a 
nucleus), Nucleus (nucleus in the neutron sea), Vortex (vortex in uniform matter), Uniform (uniform 
matter). For the configurations Pinned and Vortex we give the total energy Etot+AE (in MeV), 
corrected in order to obtain the same number of neutrons as in the configurations Nucleus and 
Uniform. We also give the cost E vor (in MeV) to create the vortex pinned on the nucleus and the 
vortex in uniform matter, and finally the pinning energy E p i nn i ng (in MeV). 
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-3.85 
















Vortex 


1.40 


117.20 


-172.05 


120.66 


124.78 


2.61 


















Uniform 


1.35 


120.15 


-198.97 


118.33 















SkM* 


noo=0.004 fm" 3 








A 


N 


Epair 


Etot 


Etot+^E 


£jvor 


-Spinning 


Pinned 


2.05 


378.80 


-413.59 


-626.93 


-499.45 


13.70 


















Nucleus 


2.00 


440.98 


-524.10 


-513.16 






-1.63 
















Vortex 


2.05 


278.08 


-440.18 


417.50 


429.03 


15.33 


















Uniform 


2.00 


283.70 


-495.26 


413.70 
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SkM* 7100=0.008 


fm" 3 

1111 








\ 

A 


N 


TP 


TP 

^tot 


77, , _i_ A 7? 


rpvor 


TP ■ 

Spinning 


_r innecl. 


o.UO 


( oO.Uo 


-919.51 


459.09 


oiz.oy 


39.67 


















Nucleus 


3.00 


802.57 


-1085.98 


472.92 






3.95 
















Vortex 


3.05 


659.10 


-1002.84 


1459.46 


1486.87 


35.72 


















Uniform 


3.00 


668.09 


-1106.10 


1451.15 















SkM* n 


00=0.012 fm" 3 








A 


N 




Etot 


Etot+^E 




E ■ ■ 

Spinning 


Pinned 


3.62 


1057.68 


-1235.92 


1420.66 


1436.74 


57.07 


















Nucleus 


3.53 


1062.12 


-1412.41 


1379.67 






6.38 
















Vortex 


3.60 


926.77 


-1317.69 


2401.44 


2398.33 


50.69 


















Uniform 


3.53 


925.91 


-1433.06 


2347.63 









SkM* noo=0.016 fm" 3 




A 


N 




Etot 


E tot +AE 




-Spinning 


Pinned 


4.25 


1397.33 


-1534.87 


2827.25 


2970.22 


73.78 


















Nucleus 


4.20 


1430.97 


-1749.45 


2896.34 






7.97 
















Vortex 


4.25 


1280.31 


-1636.09 


3864.68 


3901.98 


65.81 


















Uniform 


4.20 


1289.09 


-1789.18 


3837.70 
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SkM* 7100=0.025 fm- 3 




A 


TV 




TP 

mot 


TP _1_ A TP 


rpvor 

Hi 


Spinning 


Pinned 


5.65 


2261.44 


-1897.82 


7286.06 


7358.10 


116.90 


















Nucleus 


5.60 


2274.19 


-2161.68 


7241.19 






16.33 
















Vortex 


5.65 


2128.52 


-2026.77 


8249.63 


8286.71 


100.57 


















Uniform 


5.60 


2135.08 


-2227.12 


8186.13 















SkM* 


noo=0.038 fm" 3 








A 


N 


E ■ 

■'-'pair 


Etot 


Etot+^E 


■npvor 


E ■ ■ 

Spinning 


Pinned 


7.55 


3481.75 


-1857.15 


15563.57 


15566.52 


154.43 


















Nucleus 


7.5 


3482.14 


-2130.98 


15412.09 






18.27 
















Vortex 


7.55 


3342.92 


-1965.67 


16476.96 


16486.52 


136.16 


















Uniform 


7.5 


3344.18 


-2210.59 


16350.36 









Table IV: The same as in Table EH but for the SkM* interaction. 
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V. COMPARISON WITH PREVIOUS MODELS 



The density dependence of the pinning energy presented above is remarkably different 
from that obtained in previous works. 

Epstein and Baym [10] calculated the pinning energy solving the Ginzburg-Landau equa- 
tion. They assumed that in the pinned configuration 'the order parameter in the nucleus 
would be similar to that for the unperturbed vortex line'. They calculated the kinetic energy 
associated with the vortex axis at a distance s from the center of the nucleus, assuming the 
same functional form and spatial dependence for the order parameter associated with the 
vortex, with or without the presence of the nucleus. Due to their basic assumption, the 
difference in energy between the pinned (s = 0) and the interstitial (s — ► oo) configurations 
only involved the condensation energy (cf. Eqs. (5.2), (A13) and (C3) of their paper). As 
a consequence, the pinning configuration becomes favoured when less pairing energy is sup- 
pressed by the vortex, than in uniform matter. This in turn implies, that the pairing gap 
inside the nucleus must be smaller than in the outer gas: this happens in their calculation 
for all densities larger than 10 12 98 g cm -3 , or 0.0055 fm -3 . 

The results obtained by Epstein and Baym are qualitatively similar to those obtained 
in the semiclassical model studied by Donati and Pizzochero using the Argonne force, 
who found antipinning at low densities and pinning for densities larger than about 0.01 
fm -3 (cf. Fig. [281(c) and Fig. [33]). As we discuss in Appendix C, the differences between 
the quantal and the semiclassical results are not due to differences in the definition of the 
pinning energy. In order to have some better insight, we recall that the pinning energy 
is the difference between the excitation energy associated with a vortex created in 

uniform matter and the excitation energy E^L associated with a vortex built on the nuclear 
volume (cf. Eq. flSJ)). In the model of ref. [11] these excitation energies result from two 
essential contributions. One of them is the condensation energy E condl associated with the 
pairing fields calculated with and without the presence of the vortex. The other is the 
kinetic energy Ej\ ow associated with the velocity field created by the vortex. For each of the 
four relevant configuration (pinned, nucleus, vortex, uniform), Ej\ ow and E cond are obtained 
integrating respectively the kinetic energy density 1/2 $ 2 /n, and the condensation energy 
density -3/8 [n(p, z)A 2 (p, z)/(X — U HF (p, z))]. The excitation energy associated with the 
pinned vortex is estimated as E%£ c = [E flow + E cond ] pinned - [E flow + E cond ] nudeus . Similarly, 
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E uZf = \ E fiow + E cond ] vortex - [E f i ow + E cond ] uniform . Finally, the pinning energy is obtained 
as E pinning = E™ r c - E™ r if . The contribution of E fiow to the pinning energy is then E fiow = 
[Efiow]pinned — [E fi ow ] nucleus — [E fi ow \ vortex + [E f i ow ] uniform, and similarly for E cond . The values 
of Efiow and E cond obtained in the semiclassical model of Ref.fll] are shown in Fig. [28](c) as 
a function of density. It is seen that E cond is negative for densities larger than about 0.01 fm 
~ 3 . This is associated with the fact that the pairing gap in the nucleus is smaller than in the 
external neutron gas. Consequently, creating a vortex passing through the nucleus destroys 
less pairing energy than creating a vortex in uniform matter. The gain in condensation 
energy is largest, when the pairing gap in the neutron gas is maximum {kp ~ 0.8 fm -1 , or 
n ~ 0.02 fm -3 , cf. Fig. [3]). The radius of the vortex core is 2-3 fm larger in the pinned than 
in the uniform case (cf. Fig. [29]); this leads to a modest reduction of the kinetic energy in 
the pinned case. As a consequence the contribution of Ef\ ow is negative and rather small. 

In the semiclassical model, the radius of the vortex core is determined requiring that 
the absolute value of the condensation energy density be equal to the value of the kinetic 
energy density associated with the vortex flow. In uniform matter the vortex core is equal 
to the coherence length £ which is of the order, or larger than the nuclear radius Rq . At 
high neutron densities in the crust £ > R ; the nucleus is included in the vortex core and 
cannot influence its properties. The semiclassical model is then very similar to Epstein and 
Baym's and the pinning energy is only due to the contribution from the condensation energy 
(cf. Fig. 1281(c)). At lower neutron density, £ ~ R . In this case the kinetic energy density 
increases on the nuclear surface leading to an increase of the core radius of 2-3 fm. In this 
case the pinning energy receives a contribution also from the kinetic energy. 

In Fig. [281(a) and (b) we show the quantities Ef\ ow and E cond obtained from our calculated 
potentials, pairing gaps, currents and densities associated with the SLy4 and SkM* mean 
fields. For densities lower than about 0.02 fm -3 , the results obtained with the two mean 
fields are rather similar, and are substantially different from the semiclassical results: both 
Efi ow and E cond are much larger in absolute value and E cond takes positive values. This 
difference is directly related to the much larger vortex radius we calculated in the pinned 
configuration (cf. Fig]25l). In turn, this is associated with shell effects, which have been 
discussed above in Sec. IV and which are not considered in the semiclassical model, As a 
consequence, in our calculations pairing is suppressed not only inside the nuclear volume, 
but also in a large layer (6-7 fm thick) in the interface region between the nuclear surface 
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and the neutron gas (cf. Fig. [15]). This changes the energy balance as compared to the 
semiclassical model, and makes it less favourable to build the vortex on the nucleus. One 
should also add that, due to proximity effects, the pairing field does not vanish in the nuclear 
volume, as it happens in the Local Density Approximation with the Argonne potential (cf. 
Fig. IHliTOl) J27]. On the other hand, the large core radius implies the suppression of the 
velocity field in a large volume (cf. Fig. [16]), leading to large negative values of Efi ow . 

At densities larger than about 0.02 fm -3 , one can notice a change in the density depen- 
dence of Efi ow and E con d in the case of the SkM* interaction: Efi ow increases and E con d 
decreases, reaching values closer to the semiclassical ones. This is related to the weakening 
of shell effects for A > 10 MeV with this interaction (cf. Fig. 1221) . 

While in the semiclassical model the pinning energy is essentially the sum of Ej\ ow and 
E con ^ in our self-consistent calculations, the changes in the mean field induced by the vortex 
play an important role. This can be seen in Figs. [281(a). (b). where we compare the pinning 
energy with the sum Ej\ ow + E con d. The contributions of the changes in the mean field to 
the pinning energy are of the order of 5-10 MeV: they are smaller than the range of values of 
Efiow o r E con d but are relevant for a quantitative assessment of the pinning energy, leading 
in some cases to a change of its sign. 



VI. CONCLUSIONS 



Vortices are one of the clearest examples of the quantal nature of superfluids. In particular 
of Fermi superfluids made out of a large number of overlapping Cooper pairs. Rotation can 
only happen by giving one [y = 1) or more [y = 2,...) units of angular momentum to 
each Cooper pair. It is then not surprising that the properties of these rotonic excitations 
in a rather complex many-body system like that of finite nuclei forming a Coulomb lattice 
immersed in a sea of free neutrons display a complex dynamics, reflecting simultaneously the 
variety of quantal behaviour associated with both finite and infinite systems. In particular 
spatial quantization typical of finite many-body systems which provides special bunchiness to 
the single-particle levels and is responsible not only for the special stability of certain nuclear 
species, but also determines the structure of single-particle resonances in the continuum: 
nuclear structure and nuclear reaction aspects of the same phenomenon. 

The effective mass (m k ) associated with the nucleon-nucleon interactions used to describe 
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Figure 28: Pinning energy, flow energy and condensation energy (in MeV) derived from the HFB 
calculations using the SLy4 mean field (a) or the SkM* mean field (b). In (c) we show the same 
quantities calculated in the semiclassical model. 



40 




Figure 29: The radii of the free and pinned vortex according to our quant al model (cf. Fig. [25]) 
calculated with the SLy4 and with the SkM* mean fields are compared with the radii obtained in 
the semiclassical model of ref. [HI] 

the system makes these effects more or less important depending on whether m& is smaller 
than, or close to the bare mass m, in keeping with the fact that the associated mean fields 
are deeper or shallower, leading to larger or smaller energy separations between resonances 
of different parities, respectively. 

While the associated results for the pinning energy are quite different, in the first case 
(m k < m) producing a bell shape like pinning energy behaviour as a function of density, 
and in the second case (m^ ~ m) predicting a linearly increasing behaviour, both are quite 
different from the semiclassical result. Consequently, while we cannot say which of the two 
behaviours is correct due to the existing latitude in the parametrization of effective nuclear 
forces, one can conclude that there is likely no alternative to a full quantal treatment of 
the nucleus-vortex interaction in the determination of the pinning energies. To become 
quantitative, one needs to consider on par with the bare nucleon-nucleon interaction, the 
induced interaction arising from the exchange of density and spin modes, taking into account 
both self-energy, vertex and induced interaction type of processes. There is a formidable task 
lying ahead, if one wants to become quantitative and realistic concerning the vortex-nucleus 
interaction in the situation typical of the inner crust of a neutron star. We believe that one 
can attack this problem with good chances of success capitalizing on the results obtained in 
the present work. 
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VII. APPENDIX A 



The HFB equations can be written in cylindrical coordinates as 

(T + U HF (p, z) - \)Uqrn(p, Z, (j)) + A(p, Z, (f))v qm (p, Z, 0) = E qm U qm (p, Z, 



(13) 

A*(p, 0, z)u qm (p, z, (f>) — (f + U HF {p ) z) - \)v qm (p, z, (/)) = E qm v qm (p, z, (/)). 



The kinetic energy operator, f = -yV ^^ V contains the effective mass p(p, z) associated 
with the adopted Skyrme interaction, while A is the pairing field and U HF is the mean 
field potential calculated using the Skyrme functional. We include all the terms of the 
Skyrme functional listed in Eqs.(2.1) of ref. |23|, except for the spin-orbit and the spin 
gradient terms. The index q labels the different quasiparticle energies E and amplitudes u ) v 
associated with a projection m of the orbital angular momentum along the z— axis. Since 
we neglect the spin-orbit interaction, the spin degree of freedom is taken into account simply 
by the degeneracy factor g = 2. 

For a system with axial symmetry respect to the z-axiz, the characteristic ansatz for the 
study of the vortex is j^jj, 25] 

A(p,z,<f>) = A(p,z)e^, (14) 

where v — 0, 1, 2... is the vortex index. 

The normal and abnormal densities are given by 



7l(p, Z, (j)) = 2 V qrn(p, Z, 0>* m (p, Z, (j)) (15) 

qm 

and by 

K(p,z,</>) = ^u qm (p,z^)v* qrn (p,z^), (16) 

qm 

and the gap is given by 

A(p, z, (j)) = -V pair K(p, z, (j)). (17) 

The Hartree-Fock field U HF is calculated according to Eq. (2.8) of ref. |23||. We assume the 
following form for the dependence of the quasiparticle amplitudes on the angle 0: 

(18) 

v qm (p,z,(f))~e i( - m -^. 
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In this way, Eqs.(fT6l) and ( fT71) lead to a pairing field which is consistent with the ansatz 

We shall only solve the axially symmetric equation for the neutrons. According to the 
study of Negele and Vautherin [l|], who neglected pairing in the minimization of their en- 
ergy functional, the favoured number of protons is mostly determined by shell effects, and 
varies from Z = 40 to Z = 50. Because we neglect the spin-orbit interaction, the natural 
choice for us is to assume Z = 40, corresponding to a shell closure. However, more recent 
studies including pairing have shown that the value of Z can depend on the details of the 
functional, being sometimes rather different from Negele and Vautherin's solution [18||. We 
have not investigated whether a different choice of Z would lead to important changes in 
our results. The protons are deeply bound and for this reason they are not expected to 
be influenced significantly by the vortex, and thus the proton field should essentially keep 
spherical symmetry. We then perform a simple spherical HF calculation for the protons, 
enclosing them in a spherical box of radius 15 fm, and carry out a spherical average of the 
neutron part of the potential felt by the protons (in particular at each iteration we replace 
the neutron density n(p, z) by an isotropic n av (r) obtained averaging n(p, z) on the surface 
of a sphere of radius r). 

We shall solve the HFB equations in a cylindrical cell of radius pb ox and height h^ ox (cf. 
Fig. BJ, imposing the boundary condition u qm {p^ z ) (j)) = v qm (p,z,cj)) = for p = pb ox , or 
z = ±hb ox /2. We shall expand the quasiparticle amplitudes on a single-particle basis writing 

U qm (p, Z,<t>) = ^ U hi;qm 0hm(p) 4>i{z) (19) 
ih 

i(m—v)cf) 



Vqm(p, Z,</>) = Vhi-gm (j) h{j n-v){p) ^i(z) ^ 



ih 

In the z— direction we have taken the set of standing waves 



/ 2 7T 2tt 

^i(z) = \ 7 — sin(fci(z + h box /2)), k { = - — , - — , ... (20) 

V >^box i^box >^box 

In the p direction we have taken the set of solutions associated with the Schrodinger equation 
for free neutrons: 



h 2 1 d d h 2 m 21 



-P-T + 



km(p) €hm<l>hm(p), (21) 



2p p dp r dp ' 2p p 2 

where /x indicates the bare nucleon mass. Projecting the HFB equations on the cho- 
sen single-particle basis, we obtain two algebraic equations for the expansion coefficients 
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^2hi(^ lhi j^ m ~ tyUhi;qm + Y2hi ^lhij,m^hi;qm + Y2hi ^Ihij ,mVhi;qm ~ E qm Uij ;qm (22) 



^2ihi ^hlij,mUhi;qm S/ii tyVfii,qm ^hi^lhij,m-u^hi;qm ~ E qm Vij ;qm . 

The matrix elements of the kinetic energy, of the mean field and of the pairing field are given 
by 



K lhijjm = 27rJ dpdzp (j)i m (p)ipj(z) 



(j>hm(p)^i(z)] 

(23) 



p(p,z) 2p(p,z)J 2 p(p,z) 
Um hm = 2ttJ dpdzp (j> lm (p)^(z)U HF (p, z)0 hm ^(z); (24) 

A^y jm = 27r J rfprfzp 0/ m (p)^j(^)A(p,z)^ (m _^)^(z). (25) 

The total energy is given by E tot = E HF + E pair , where is E HF calculated according to 
Eq. (2.2) of ref. |23| and E pair is obtained by 

E P air = ~\ J dpd^ pK*(p, z, 0)A(p, z, (f)) (26) 

In the literature, one finds different choices concerning the boundary conditions applied 
to calculations in Wigner-Seitz cells. In Fig. [30] we compare the values of the gap at the 
Fermi energy Ap obtained with the SLy4 mean field as a function of the Fermi momentum in 
a spherical box of radius R= 30 fm, with those independently calculated in infinite neutron 
matter and already shown in Fig. [31 We have used three different boundary conditions for 
the single-particle wavefunctions: (a) all the wavefunctions vanish at the edge of the cell 
(the boundary condition adopted in this paper) (b) the wavefunctions of even orbital angular 
momentum and the derivatives of the wavefunctions of odd orbital angular momentum vanish 
at the edge of the cell (the boundary condition adopted by Negele and Vautherin); (c) the 
same as (b), but exchanging the role of odd and even angular momenta. It is seen that 
the agreement with the result obtained in infinite matter is of the same quality in the three 
cases (cf. also the discussion in ref. [28|]). 

We have also checked that with our boundary conditions (a) we reproduce accurately the 
values of the energy per particle in neutron matter calculated analytically with the Skyrme 
interactions we have considered. As we discussed in Section IVB, using boundary conditions 
(a) we have also been able to reproduce the pairing gap associated with the vortex (v — 1), 
calculated by Bulgac and Yu in infinite neutron matter. 
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For calculations with protons, that is, with a nucleus at the center of the box, we have 
found that in all the calculations reported in this paper, the pairing gap at large distances 
(p > 20 fm) resumes the value calculated in uniform neutron matter at the relevant density. 
Obviously with boundary conditions (a) the pairing gap drop quickly to zero at the edge 
of the box. We have checked that this does not influence the calculations of the pinning 
energy, which is obtained subtracting the energies of different configurations, as discussed 
in the next Section. 

In all our calculations we have used boxes with radii equal to, or larger than, 30 fm. In 
other works, where one has to calculate the pairing gap inside the actual Wigner-Seitz cells, 
one finds some clear dependence on the boundary condition for cells smaller than about 20 
fm jisj]. In this case, one should consider in detail the influence of the Coulomb lattice 

3 
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Figure 30: Comparison between the gap obtained requiring that the radial single-particle wave- 
functions ul(t) vanish at the edge of the spherical box for all orbital angular momenta L (crosses), 
or that ul{R) (for L even) and duL/dR{L odd) vanish at the edge of the box (stars), or that ul{R) 
(for L odd) and duL/dR{L even) vanish at the edge of the box (circles). The solid curve shows the 
pairing gap calculated in infinite neutron matter. 

VIII. APPENDIX B 

In this Appendix we give some details concerning the numerical accuracy of the calcula- 
tion of the pinning energy. This is important, because the pinning energy, which is of the 
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order of a few MeV, results from the difference of the total energies of the various configura- 
tions, which are of the order of a few GeV. The main parameter that influence the accuracy 
of our calculation is the mesh size. Furthermore, we have to check that the size of the box 
( that is, its height hb ox and its radius pb ox ) is sufficiently large, to obtain stable results. 
As we discussed in Section IVC, we expect that values like hb ox = 40 fm and p^ ox = 30 fm 
are sufficiently large to guarantee that the pinning process is not significantly influenced by 
the boundary. In fact, the typical vortex radius, as measured by i?go%, is about 15 fm (cf. 
Fig. [25]). Moreover, the calculation of the pinning energy is more sensitive to the value of 
Pb ox than to h\) OX ^ because the distortion of the pinned vortex is much more relevant on the 
equator (cf. for example Fig. [14]). We have studied in a few cases the dependence of the 
computed values of the pinning energy on p5 OX , making calculations of the relevant configu- 
rations with different boxes. In Tab. [V] and in Tab. [VH we show the values of the pinning 
energy calculated with different box sizes for the SII and SLy4 interactions. The values 
show very modest fluctuations with the box size for the points associated with densities up 
to about n < 0.017 fm -3 , for which we can estimate an uncertainty of the order of 200 keV. 
For the highest density points, the fluctuations become larger, but do not exceed 1-2 MeV. 
This is related to the fact that the vortex radius increases with density, particularly with 
the SII and SLy4 interactions (cf. Fig. 1251) . 



A [MeV] 


rioo [fm -3 ] 


Pbox= 30 fm 


p box = 25 fm 


3 


0.004 


-1.51 


-1.65 


5.8 


0.011 


3.26 


3.46 


7 


0.016 


1.64 


1.34 


11.3 


0.035 


2.34 


1.77 



Table V: Pinning energy (in MeV) calculated with the SII interaction at various values of the 
chemical potential and of the asymptotic neutron densities, for two different radii of the cylindrical 
box. 

The size of the mesh in p and z used in the code that solves the HFB equations and 
calculates the energy is crucial for the accuracy of the calculations. We have studied the 
influence of the mesh (for the p coordinate) in the case of the SLy4 interaction at the density 
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A [MeV] 


rioo [fm~ 3 ] 


Pbox= 35 fm 


Pbox= 30 fm 


Pbox= 25 fm 


5.7 


0.011 


0.45 


0.42 


-0.3 


7 


0.017 


6.07 


6.21 


5.54 


10 


0.037 


1.64 


-0.41 


1.78 



Table VI: Pinning energy (in MeV) calculated with the SLy4 interaction at various values of the 
chemical potential and of the asymptotic neutron densities, for three different radii of the cylindrical 
box. 
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Figure 31: Pinning energy calculated with the SLy4 interaction at the density = 0.001 fm -3 as 
a function of the mesh size used in the p-coordinate. 

rioQ = 0.0013 fm -3 and the results are shown in Fig. [3U The fluctuations of the pinning 
energy for a mesh smaller than 0.3 MeV are of the order of hundreds of keV (< 0.2 MeV). 



IX. APPENDIX C 

In a semiclassical study of the problem under discussion jll] the pinning energy was 
defined as the difference of the energies associated with the pinned and the interstitial 
configurations, taking explicitly into account the presence of one of the neighbouring nuclei 
(cf. the left part of Fig. [32). This definition of pinning energy reduces to ours, if the 
interaction of the vortex with the second nucleus is negligible. We have tested the importance 
of the neighboring nucleus within the semiclassical model, for the range of densities we have 
considered. We first computed the energy of each of the relevant configurations using the 
semiclassical theory, the two pairing interactions (Argonne and Gogny) and the parameters 
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Figure 32: (left) The scheme used by Donati and Pizzochero [11] where the pinning energy takes 
into account the presence of a neighboring nucleus; (right) the scheme we used, where the pinning 
energy is evaluated as the difference between the costs to excite a vortex on a nucleus and in uniform 
matter. 



presented in Ref. [11] , instead of solving the HFB equations. We then calculated the pinning 
energy using the scheme proposed in Ref. [11]: the resulting values are connected by the 
dashed curves in Fig. [33], and they reproduce quite accurately the values reported in ref. fll|. 
already shown in Fig. [281(c) in the case of the Argonne interaction. In Fig. [33] we also show 
the pinning energies calculated using the semiclassical energies, but using our scheme (cf. 
the right part of Fig. [32]). The two schemes yield an almost equal result, showing that the 
differences existing between the quantal and the semiclassical calculations are not due to 
the geometrical scheme used to derive the pinning energy, and that the contribution of the 
second nucleus plays in fact a marginal role, also following the formalism of ref. [11]. The 
situation may change at larger densities, where the vortex radius becomes of the order of 
the internuclear distance, as we have discussed in Section IV C. 
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Figure 33: Solid and dashed curves show respective! 
scheme (Fig. right part), or the scheme of Ref. 
have computed the energies of the different configurations using the semiclassical model, with the 
Argonne pairing interaction (left panel) or the Gogny pairing interaction (right panel). 
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